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We study with lattice techniques the locahsation of gauge fields on domain wall defects 
in 2+1 dimensions, following a scenario originally proposed by Dvali and Shifman for 3+1 
dimensions, based on confining dynamics in the bulk. We find that a localised gauge zero- 
5^ \ mode does exist, if the domain wall is wide enough compared with the confinement scale 

in the bulk. The range of applicability of the corresponding low-energy effective theory is 
determined by the mass gap to the higher modes. For a wide domain wall, this mass gap is set 
by "Kaluza-Klein modes" as determined by the width. It is pointed out that in this regime 
the dynamical energy scales generated by the interactions of the localised zero-modes are in 
fact higher than the mass gap. Therefore, at least in 2+1 dimensions, the zero-modes alone 
do not form a low-energy effective gauge theory of a traditional type. Finally, we discuss how 
the situation is expected to change in going to 3+1 dimensions. 
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1. Introduction 



In brane- world scenarios with infinite ^ |5J |^ or large ^ extra dimensions, it is assumed 
that all the fields of the Standard Model have wave functions that are localised in trans- 
verse directions, making the physics to be four-dimensional (4d) at small energy scales. The 
field-theoretical realisation of branes is related to topological defects — stable solutions of 
the classical equations of motion, which depend on the transverse coordinates only. In five 
dimensions, the necessary defect is a domain wall, in six dimensions it is a (4d) string, in 
seven dimensions a monopole, etc. (For an explicit construction of these solutions in the case 
of localisation of gravity, for instance, see Refs. 01^1 CI-) 

Small perturbations of the fields around the topological defect solutions may form a low- 
dimensional effective theory — potentially the Standard Model (for a general discussion of 
perturbations see, e.g., Ref. 0). Apart from the requirement that the wave function of the 
"fundamental" (or lowest-energy) mode be centered around the brane in transverse directions, 
such that it has long-wavelength perturbations in four directions only, a successful localisation 
poses other constraints as well. Indeed, either the perturbations of all the higher modes should 
be separated from those of the fundamental mode by a sufficient mass gap, in order not to 
be observable in low-energy experiments, or, if there is no mass gap, the perturbations of the 
higher modes should interact very weakly with those of the fundamental one. 

The possible solutions to these requirements can roughly be divided into two classes (for 
a recent review, see Ref. [5]). On one hand there are mechanisms which, in one way or the 
other, invoke effects related to gravity (such as a "warped" metric ^lOl)- In this way scalars 
[niII21, fermions [H], |I3]-|ini, Abelian gauge fields HHH, and gravity can be 

localised on a brane, although the construction of the full Standard Model is still far from 
being achieved. On the other hand, there are mechanisms which work in flat spacetime and 
are purely field theoretic in origin. For example, the localisation of fermions on a brane may 
be due to the presence of fermionic zero-modes on topological defects 012111, while scalars 
can be attached to the brane through their interactions with the field forming the topological 
defect (and the perturbations of the defect forming field itself constitute a light scalar field). 

The most difficult problem along these lines seems to be the localisation of massless non- 
Abelian gauge fields (for a general discussion, see Ref. [U|). If the mechanism related to 
gravity is used, for instance, then for typical warp factors, the bulk gauge fields have a spec- 
trum which is not separated by a mass gap from the localised modes, so that no effective field 
theory can be constructed.^ At the same time, in the case without gravity, no perturbative 
mechanism of gauge field localisation is available at present, as far as we know.^ 

A very interesting non-perturbative purely field theoretic idea for gauge field localisation 

^Both massless and massive vector fields can be localised on a brane, such that the fundamental mode is 
separated from the higher ones by a non-zero mass gap, if the warp factor is tuned appropriately "by hand" 
[2211211 (see also Ref. 

^Although of significant interest, we do not consider mechanisms related to a high degree of unbroken 
supersymmetry in this paper (see, e.g., Ref. |25| and references therein). 
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was put forward some time ago by Dvali and Shifman [^Hl- It uses the fact that non-Abehan 
gauge theories are strongly coupled in the infrared and that a mass gap can be generated 
for the bulk modes by confinement effects. In short: consider a confining theory in the bulk, 
based on some group G. Construct a topological defect that "eats up" the necessary number 
of dimensions and breaks this symmetry down to G' inside the defect, while leaves it intact 
outside. Then the gauge fields related to G' are localised on the defect and are separated by 
a mass gap from the bulk modes, which are massive because of confinement. The resulting 
low-energy effective theory is thus a four-dimensional gauge theory based on G'. 

Of course, to use this idea for the construction of a realistic theory, one would first have to 
find a confining gauge theory in 4 + n dimensions, with n > 1. Non-Abelian gauge theories 
of the usual type are, however, not renormalisable when extra dimensions are involved. In 
five dimensions, for instance, lattice simulations [^Tj do not reveal any second order phase 
transition that can be used for a continuum formulation of a confining theory.^ Thus, to 
have a simple renormalisable framework, Dvali and Shifman assumed the bulk dynamics to 
be that of 4d gauge theory, so that the low-energy dynamics is that of three-dimensional (3d) 
gauge theory. 

Even if plausible after this simplification, the idea does still involve some untested assump- 
tions. First of all, it may not be immediately clear why the non-perturbative confinement 
effects, acting in the bulk and admittedly suppressed on the brane, could not generate an 
effective mass term for the brane gauge fields; this depends after all on the precise boundary 
conditions that the bulk phase poses on the localised modes UHl- Another question, neces- 
sary for understanding whether the effective theory is truly lower dimensional than the bulk 
theory, is related to the magnitudes of the confinement induced mass gap Mc, the typical 
energy scale of the low-dimensional theory nii, and the masses of higher "Kaluza-Klein" (KK) 
excitations M^, localised on the brane. ^ To hide extra dimensions at small energies, the scale 
mi must be much smaller than the masses of particles in the bulk as well as those of the 
KK-excitations: mi <C M^, mi <C Mc- Otherwise, the effects of higher dimensions are not 
suppressed. This point has not been discussed in Refs. [5^11131) as far as we can judge. 

Given the open ends, it is the purpose of this paper to test the idea of Ref. [2^1 with simple 
non-perturbative lattice simulations. We consider SU(2) gauge theory coupled to an adjoint 
scalar field with an effective space-dependent mass parameter induced by the kink solution, 
as in the original proposal 26 . Since going to higher dimensions leads us to a shaky ground 
we will deal, following Ref. |2^, with renormalisable theories only, i.e. consider at most a 4d 
bulk. For better numerical resolution and since the physics arguments are almost unchanged 
(see below), we will however reduce the dimensionality to be three in total, such that the 
low-energy effective theory is supposed to be two-dimensional (2d). 

Our lattice simulations do support the Dvali-Shifman conjecture on the existence of a gauge 

''For recent discussions and references see, e.g., Refs. |28| . 

*Even though there is no compact dimension involved, we refer to the tower of localised states above the 
fundamental mode as Kaluza-Klein excitations, due to the fact that they have a similar spectrum. 
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field localised on the brane. In other words, if the dynamics of the theory is probed with 
external sources separated by large distances along a (wide enough) brane, a force specific to 
a massless localised gauge mode will indeed manifest itself. We find, however, that in this 
limit there is no hierarchy of (dynamical) scales between the 2d localised mode and certain 
higher excitations and, therefore, that the low-energy theory cannot be considered to be a 
traditional 2d gauge theory. 

This paper is organised as follows. In Sec. El we review the basic mechanism as proposed 
by Dvali and Shifman 26 . In Sec. |31we derive its signatures for a specific observable, the 
static force between two heavy test charges living on the brane, and in Sec. ^ discuss how 
the signatures change if the world- volume is finite. In Sec. [S] the system is put on the lattice 
and the observables are written in a form accessible to Monte Carlo simulations. The main 
results are presented in Sec. and our conclusions as well as a brief outlook, in Sec. [Tj 



2. Overview of the Dvali-Shifman mechanism 

The model suggested by Dvali and Shifman has two scalar fields, a gauge singlet r] and an 
adjoint scalar X) and has the following Euclidean action in the case of a 3d bulk: 



Se = j d^xCE , (2.1) 

£e ^ ^Tr Fii + Tr [Dk,xf + A TV [x^]^ + 7 Tr [x^] (^.^ - + ^,2) 

+ \id,rjf + \^{v'-v')\ (2.2) 

where k,l = 1,...,3, Dk = dk + igAk, Ak = AIT'', x = X^T^, Fu = (l/ig)[Dfc, A], and T'' 
are the Hermitean generators of SU(A'^c), normalised as Tr \T°'T^] = 6"'^/2. Summation over 
repeated indices is understood. The coupling constants A, 7, k and the parameters and v'^ 
are assumed to be positive, and > v?. 

The classical vacuum of the theory is at x = 0) ^ = and the perturbative spectrum 
consists of the scalar singlet with the mass = 2kv'^ , the scalar triplet with the mass 
= 7^2 and massless gauge bosons corresponding to the gauge group SU(A''c). Because of 
strong coupling in the infrared, the vector boson spectrum aquires a mass gap of the order of 
the 3d confinement scale, A ~ jj-^ other words, gluons form bound states — glueballs — 
with a mass of the order of A. Depending on the relation between A and m^, bound states 
of the triplet scalar and gluons can have masses 0{h) or 0{m^). 

Now, the model of Eq. (|2.2|) always has a kink solution, x = Oi 'ni^) — t'tanh(^^^). 
By inspecting whether the fluctuation Hamiltonian around this solution develops a negative 
eigenvalue, it is seen that the solution is unstable against x-field condensation near z = 0, 
provided that 

■yv"^ > -I- -m^niri . (2.3) 
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In this case the stable classical solution contains a non-zero x held as well, and has the 
asymptotics Xs ~^ 0) ^ ^ at z — > ±00, whereas the components xi and X2 can be chosen 
to be zero. For a general choice of parameters the explicit solution can easily be constructed 
numerically. The analytic form can be found for a specific choice of parameters, namely for 



2^n' = ^^^v^ , (2.4) 



and is given by 



ri{z) =vi&n\i{m^z) , x^iz) = J- — ^ — . (2.5) 

y 7 — A cosh(m^z) 

In the special case k = 7 = A, 'q{z) remains the same while 

= ^^511^; . (2.6) 

cosh(m^2:j 

To get the mechanism to work, the following choice of parameters is proposed. Inside the 
defect the SU(A'^c) symmetry is partially broken to G', and the masses of the corresponding 
vector bosons are assumed to be large compared with the confinement scale A, 

7(f;2 - n2) 2 



A 



» . (2.7) 



At the same time, the width of the domain wall is chosen to be larger than the inverse of 
the confinement scale, ^ m^. These two requirements suppress the non-perturbative 
confinement effects inside the defect. Together with Eq. (|2.l-i|l . they pose restrictions on the 
parameter space, but given that there are several parameters at our disposal (it, f;. A, 7, k), 
all the requirements can easily be satisfied simultaneously. 

At the quadratic level the spectrum of perturbations around the domain wall contains a few 
normalizable localised scalar modes (there is one zero- mode associated with brane translations 
and fluctuations), a continuous spectrum of scalar excitations that starts from min(m,y, m^^), 
and a gapless continuous spectrum of vector excitations corresponding to the unbroken group 
G'. Now, because the true non-perturbative spectrum of bulk gauge excitations is massive, it 
was conjectured in Ref. [^H] that the vector bosons of the unbroken group cannot escape the 
brane and thus the true spectrum of gauge excitations around the domain wall does contain 
normalisable vector modes related to G'. If true, the low-energy effective theory is just a 2d 
gauge theory in our case. 

Further support for the Dvali-Shifman idea has been provided by Arkani-Hamed and 
Schmaltz, who argued j^H, based on the 't Hooft - Mandelstam picture of conflnement, 
that the bulk acts as a Neumann boundary condition for the gauge fields corresponding to 
G', such that there is a true localised zero-mode, unlike in the case of a Dirichlet boundary 
condition, whereby a zero-mode is excluded.^ Since the physics of confinement is involved, 

^For further work on the topic see, e.g., Ref. 
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however, a non-perturbative check of the mechanism would be welcome, and this is one of the 
aims of the present work. Another problem, related to mass gaps and the dimensionality of 
the low-energy theory, has been mentioned already in the Introduction, and will be explained 
in more detail below. 



3. How to probe the properties of the low-energy effective theory? 

Since our main interest here is in the gauge fields, we will ignore in the following all the 
dynamics related to the scalar singlet field 77. To achieve this we go to the rest frame of 
the domain wall, and also treat it as infinitely rigid, which means that we remove the 3d 
zero-mode related to translations and fluctuations of the domain wall. 

The remaining dynamical degrees of freedom constitute the SU(A'c) gauge + adjoint Higgs 
theory, formally defined by the action 

Se = J d^xCE , (3.1) 
Ce = ^TrFii+Tr[Dk,xf + m^iz)Tr[x^]+X{Tr[x^])\ (3.2) 

where m'^{z) is simply some profile for the mass of which we fix "by hand" to be 

where, in the notation of Sec. |2 m2 = = represents the mass of the scalar triplet 
outside the brane, while = j{u'^ — v'^) is negative and is related to the mass of the scalar 
boson inside the brane. As already discussed (cf. Eq. H2.7() 'l. we assume that the masses of 
the vector bosons inside the brane are large, niyy ~ —g^m^/X S> g^. The domain wall width, 
^, is for generality now treated as a parameter independent of m^, and is supposed to be 
large enough, 

^ » ^ , (3.4) 

to suppress the influence of bulk confinement on the localised gauge field. In addition, to 
have condensation of x inside the brane, we must require (cf. Eq. ()2.3|) ) that m\ < —m2/i- 

The dimensionless combinations of the parameters, determining the actual dynamics, can 
be chosen as 

a = ig , x = ^ , yi,2 = , (3.5) 
9^ 9* 

where the mass parameters (the only ones requiring renormalisation in three dimensions) are 
for convenience assumed to be evaluated in the MS scheme at the scale fl = g^. 

Let us assume now that the Dvali-Shifman conjecture is correct and estimate the parameters 
of the localised gauge theory. Since the wave function of the massless mode is localised on 
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a length of the order of the width of the domain wall, the 2d effective gauge coupling g2 is 
simply 

1 ~ 4 • (3-6) 

If two opposite static test charges are put on the brane, one would then expect that they are 
attracted, at large distances (to be specified presently) with a force F ~ gr^ that does not 
depend on the distance. Note that this holds also for G' = U(l), since the Coulomb potential 
is linear in 1+1 dimensions. 

Besides the massless mode one would expect to also have a whole tower of states with 
the same quantum numbers and with an energy spacing of the order of ~ 1/^, with E 
coming simply from the uncertainty principle. These states are not seen at distances r ^ £ 
or at energies E <^ so that the force derived above should be valid for r ^ i. For 
^ ^ ^ the 2d Coulomb law is expected to be replaced by the 3d Coulomb law with 
the force F ~ while at even smaller distances the massive W will contribute as well, 

changing the numerical coefficient in front of the 1/r-dependence of the force. 

If, on the contrary, the localised zero-mode acquires a mass m-y due to interactions with 
the bulk modes, the force between the test charges at r ^ i will have a Yukawa character, 
F ^ 02 exp(— m-yr), whereas for smaller r the behaviour is still as described above. 

These two different behaviours of the force at r ^ £ can be distinguished in lattice simu- 
lations, as will be described below, so that a conclusion can be reached on the existence of a 
localised vector zero-mode. As we will see, the presence of a vector zero-mode is confirmed, 
and the dependence of the 2d coupling constant on the width of the domain wall (Eq. ()3.6() 1 
is also found to be as expected. 

Let us now discuss whether the low-energy theory can indeed be considered to be 1+1 
dimensional electrodynamics. If so, the typical energy scale Eq of the 1+1 dimensional 
theory must be smaller than the mass gap to the first excited KK mode. To be concrete, 
let us imagine adding fermions to the 3d theory and let them interact with the singlet field 
in a way that ensures the existence of fermionic zero- modes. The gauge coupling of these 
fermions is g2- Then, the theory of the zero-modes (fermionic and vector) is simply the 
Schwinger model, the solution of which tells that bosonic scalar states with the mass Eq ~ g2 
are formed. The KK tower of states decouples, provided that Eq ^ l/l. This inequality, 
together with Eq. ()3.6() and the estimate Eq ^ §2, gives i <^ g'^, which is in contradiction with 
the initial assumption of Eq. (|3.4() . In other words, we do not expect the low-energy effective 
theory to be the Schwinger model, it will rather be a more complicated 1+1 dimensional 
theory incorporating not only the massless mode but also its KK excitations. 

This expectation can also be tested on the lattice. Indeed, if for some reason the mass 
gap to the KK excitations is considerably higher than 1/i and corresponds to some distance 
scale io <^ £, then the 1+1 dimensional Coulomb law will be valid to considerably smaller 
distances. As we will see in the lattice simulations, the deviations start really from r ~ 
and are characterised by a mass gap ~ l/£. Therefore, the low-energy theory, though giving 
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an apparently 1+1 dimensional force at the distance scales 3> is in general not that of the 
zero-modes alone. 

We also note that the other logical possibility, i < does not lead to 1+1 dimensional 

electrodynamics either. In this case the spread of the potentially massless localised mode is 
given by the inverse confinement scale, in our case. Thus, the effective 2d gauge coupling 
is 92 ~ and, therefore, the typical 2d energy scale is of the same order as the mass gap to 
the bulk modes. 

To conclude this section, we note that the specific line of reasoning above is related to the 
case of 3— >2 compactification. In the last section we however provide arguments that the 
pattern is essentially the same also for 4— >3 compactification, in the generalised case that G' 
is a non-Abelian group. Unfortunately, we have nothing to say about (the most interesting) 
higher dimensional case. 



4. Finite-size scaling of the Abelian static force in a periodic box 

As explained in the previous section, the probe to be used in order to test the conjecture, 
is the static force between infinitely heavy test charges living on the brane. The motivation 
for this choice is that theoretical predictions, as reviewed in the previous section, are unam- 
biguous, and that lattice measurements, as discussed in Sec. 13 can be made rather precise, 
employing recent technical advances jHJ. To perform the lattice simulations, periodic bound- 
ary conditions will be used, so we shall first discuss how this changes the expectations for the 
static force presented above for the infinite-volume case. 

To fix the notation, let us consider the three-dimensional volume to be a Euclidean hyper- 
torus (or box with periodic boundary conditions in all directions), with a coordinate r in the 
spatial direction along the brane, t in the temporal direction along the brane, and z along 
the "bulk" , perpendicular to the brane. The extents of each direction (if finite) are denoted 
by Lr,Lt,Lz, respectively. The brane is located at z = 0. 

The static force is defined in the usual way. We introduce a rectangular Wilson loop 
W{R,T; z) in the (r, i)-plane, at some fixed z, of size Rx T. The force F{R;z) is defined 
through the potential V{R; z) as 

F(fi;.) = -li^ =^,to i-il„,y(fi.r;.) . (41) 

where we assumed that Lj = oo. 

If we imagine for a moment that the system is homogeneous and perturbative, a leading 
order computation in gives 
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where Cp is the quadratic Casimir of the fundamental representation, Cp = {N^ — l)/2Nc, 
and Pz = 27rn/Lz, with n an integer. This result is exact for an Abelian theory (Coulomb 
law), but in a non- Abelian confining theory it is only valid at short distances and, in general, 
the periodicity it displays will be lost at large distances. Nevertheless, as we will see, various 
interesting limits can be obtained from this simple expression. In the following we set iVc = 2 
as in the actual simulations, so that Cp = 3/4. 
Let us consider three limiting cases: 

Infinitely thin brane: the symmetric phase. In the confining symmetric phase, the 
force is analytically computable only at small distances (the Coulomb part). Because of the 
mass gap, we can set — > cx) in Eq. 1)4.2(1 . so that the sum J2p^ becomes an integral. Then, 

^^^'^^ = ' R«{9'CAr\ (4.3) 

= c^ig^CA? , (s'Ca)-^ « « L,/2 , (4.4) 

where a = ci{g'^CA)^ is the string tension of 3d SU(2) Yang- Mills theory with adjoint matter, 
and Ca = Nc is the quadratic Casimir of the adjoint representation. Without adjoint matter, 
the constant ci is determined numerically to be ~ 0.0281(3) in the continuum limit for 
Nc = 2 [22], and the inclusion of adjoint matter leaves the value practically unchanged j3Sl- 

Infinitely thick brane: the broken symmetry phase. In the opposite limit of a 3d 
broken symmetry phase, the three isospin components of the vector bosons convert into two 
massive vector bosons, of mass mw = gix^)^ and one massless vector boson, the photon. In 
3d, however, the photon becomes massive jHl] via interactions with monopoles [35l l36]. 

~ g'^^^rnKf" exp ) • (4.5) 

Correspondingly, the potential is computable in weak coupling only up to distances ~ m^^. 
Assuming again, for simplicity, that — > oo, we obtain 

= 4^ , m^i«i?«m-i, (4.7) 
= C2g'^m^ , <^ R<^ Lr/2 , (4.8) 

where C2 is a constant. If Lr/2 <^ Tn~^, as is realistically the case, then at large distances 
we rather encounter the fully perturbative behaviour for the photon following from Eq. (|4.2|) 
(still in the limit Lz oo), 

^2 
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Brane of a finite width. Consider finally the case of our actual interest, a domain wall of 
some finite effective width, ~ ^brane- The behaviour should now interpolate between the two 
limits encountered above. The input parameter determining the width is a = Ig^ as defined 
in Eq. (|3.5|) . Independent of a, at very small distances we still have 

^(^;^) = ^^, R-^m^. (4.10) 

What happens at large distances, on the other hand, depends on whether there is a zero-mode 
or not. If the zero-mode exists, we may expect that its contribution is according to Eq. ()4.2() . 
where is now finite and replaced with Lbrane; and we take = 0, and only one of the 
isospin components contributes so that Ci? = 3/4^1/4, 

F(°)(i2;l^l< W)~f . (4.11) 

This would be the behaviour if the confining phase outside the brane acted effectively as a 
Neumann (derivative of field vanishes), rather than Dirichlet (field itself vanishes), boundary 
condition. The first massive mode, on the other hand, contributes as 

p(i)^o.|^l<r N _ 1 sinh[mi(L,./2-i?)] 

F (i?,|.|<Lb_) - T2Lb_ sinh(m,L./2) ^^'^^^ 
- ^!^£(z!!li^, ^,L.»l, (4.13) 

4 Zi^brane 

where the mass mi is assumed non-zero. This would also be the full behaviour with Dirichlet 
boundary conditions, in which case a zero-mode is excluded. 

Let us emphasize that the linear in R behaviour of the force in Eq. 1)4. 11(1 is a characteristic 
of the Abelian theory. If the 2d low-energy effective theory were non-Abelian, the long- 
distance force would still be constant as in Eq. (|4.4j) . with exponentially small finite- volume 
corrections. This difference is caused by the finite periodic extent of the box in the R- 
direction — in infinite volume both the 2d Abelian and non-Abelian gauge theories display 
linear confinement (constant force) a la Eq. (|4.4|) . 



5. Lattice formulation 
5.1. Discretised action 

In order to test the behaviour at large distances, so as for instance to distinguish between 
Eqs. (|4.1H) and (|4.13|) . we study the system on the lattice. The discretised Lagrangian 
corresponding to Eq. (|3.2|) is 

Aatt = -j-^ y] Tr [1 - Pki{x)] 
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[Trx^la-O - T^^c x{x)Uk{x)x{x + aek)ul{x) 



+m^,,,(z)Tr[x^]+A(Tr[x^])^ 



(5.1) 



where a is the lattice spacing, Uk{x) = exp[iagAi:{x)], is a unit vector in direction k, and 
Pki is the plaquette: 



Pkiix) = Uk{x)Ui{x + aek)Ul{x + aei)u}{x) 



(5.2) 



The extents of the box are denoted by = aNr,Lt = aNt, Lz = aN^, and the lattice volume 
hy V = NrNfNz- The lattice spacing is expressed through the dimensionless combination 



2iVe 



(5.3) 



In order for physics to be the same as with the continuum Lagrangian in Eq. ()3.2|) . the 
lattice theory needs to be renormalised. In three dimensions, the only parameter including 
divergences is m?, and the divergences can be computed exactly, close to the continuum 
limit ISZIEHI: 



"^barc( 



(5m^ (//) 



- 2g^CA + X{dA + 2 
1 



(5.4) 



47ra 



+ 



167r2 



2X{dA + 2) (A - q^Ca 
4 



In— + C] - 2g'CAX{dA + 2) 



3C1 



vrS - 4(5 + /o) + 2ki - ka 



(5.5) 



where the constants C,, 5, p, ki, K4 and S have been defined in Refs. [371 EH]) C'a = -^o 
= N"^ — 1, and fi is the scale parameter of the MS scheme. The couplings g^, A require no 
renormalisation, but the approach to the continuum limit could be improved by computing 
corrections of order 0{a) 39 . As already mentioned, we choose = 5^ in the following. 
Various condensates, such as Tr[x^], also require additive renormalisation (multiplicative 
wave function renormalisation effects are 0{a) [55 ): 



dACA 
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5.2. Discretised observables 

The basic object we employ on the lattice is a Polyakov loop, 



,6 , S2 
ln— + C + — -6 
a/x 4 



(5.6) 



P{R-z) = Tt U^^I^Ut{R,nta,z 



(5.7) 
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and the correlation function of two Polyakov loops is 



C{R\ z)= (^P\nra; z)P{nra + R; z)J> . 



(5.8) 



nr=0 



The general structure of this correlator is (^31' and references therein) 



C{R; z) = wie 



LtV(R;z) 



+ W2e 



LtlV(R;z)+AE] 



+ ... 



(5.9) 



where Wi are numerical coefficients and AE is a mass gap. Therefore the static force, to be 
denoted by Fp{R;z) when extracted from the Polyakov loop correlator, can be defined as 



and the difference between Fp{R; z) and F{R; z) is exponentially small if AE > and Lt is 
large enough. On the lattice we use a discretised version of Eq. ()5.1U() . 



The reason for using Polyakov loop correlators instead of the more common Wilson loop to 
obtain the static force, is simply that a sufficient numerical accuracy is easier to reach, thanks 
to the advanced numerical techniques developed by Liischer and Weisz (and optimised 
for the present system by us). The measurement method will be described below. This 
advantage does not come without a price, however: if the system is not confining, then the 
error made by using Eq. 1)5. 1U() rather than the full expression in Eq. ()5.9() may in general not 
be small, and the correct result is only obtained after an extrapolation to the limit of a large 
Lj. Nevertheless, even data obtained at finite Lt will correctly show whether the behaviour 
of the force is constant,^ linear, or exponential in i?, thus allowing to distinguish between 



In the following, we will always refer to the object defined by Eqs. (|5.1U|) . ()5.11|) as the 
static force, while keeping the reservations just spelled out in mind, and therefore checking 
in the end explicitly for the stability of our results with respect to variations of Lt. 

5.3. Parameters and systematics 

The model of Eq. (|5.1|) has previously been studied with lattice simulations in Refs. [411 142 j . 
for the case of a homogeneous mass parameter, corresponding to a = or a = cxd. The system 
was observed to have a non-trivial phase diagram in the space of the continuum parameters 

®At large distances the behaviour of Fp deviates from a constant, due to the periodic boundary condi- 
tions. Based on the effective string picture, the behaviour expected is the same as in 2d |4()) . Fp{R;z) = 
atainh[aLt{^ — R)] for large Lr, which approaches a step function for Lt — » oo. In practice, however, the 
error bars are reasonably small only for R <^ Lr/2, where Fp can be well approximated with a constant. 




(5.10) 





Eqs. Km . Km . 
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defined in Eq. ()3.5|) : for x < Xc ~ 0.3, there is a first order phase transition separating 
the broken symmetry phase at y < i/c from the symmetric phase at y > ye- For x > Xc, on 
the other hand, there is no phase transition, meaning that the two phases are analytically 
connected. This is possible since there are monopoles in the broken symmetry phase |35[ I36j 
which replace a possible order parameter, the mass of a photon related to an unbroken U(l) 
symmetry, by a small but non- vanishing mass (given in Eq. (|4.5p ) related to a pseudoscalar 
particle |31I. The situation would be different with the gauge group SU(3), for instance, in 
which case there is an additional discrete global symmetry x ~^ ~X which gets broken at 
y = yc', thus there always is a genuine phase transition of some kind |43j . 

In terms of the parameters of Eq. 1)3. 5(1 . we choose here to study the value x = 0.20, which 
implies a first order transition at yc ~ 0.147 [32 (provided one is close to the continuum 
limit, i.e., /? ^ 1). The motivation is to have a clear distinction between the two phases, 
and to avoid complications owing to the presence of monopoles, which are more important 
in the crossover region. The parameters yi, y2 of Eq. (|3.5|) are chosen on the two sides of the 
transition, 

yi = -0.667 , y2 = 2.0 . (5.12) 
The width of the domain wall is varied in small intervals in the range a = 

V = 0.0... 5.3, 

and we also have data for a = oo. 

In order to control finite lattice spacing and finite volume artifacts, we have carried out 
simulations at two values of (3, (3 = 6,9, and with a series of volumes in the range V = 
NrNtNz = 24^... 48^. Experiences from previous studies of the same system |411 1421 OH] as 
well as from glueball computations with pure SU(2) gauge theory in three dimensions j32j 
suggest that such values are already safely in the scaling region. 

The update algorithm used for the simulations is a combination of one heat bath update 
cycle followed by four overrelaxation updates for both the gauge and the adjoint scalar fields. 
More details concerning the implementation can be found in Ref. j42j . 

5.4. Measuring the Polyakov loop correlation functions 

In a confining system at large Lt the Polyakov loop correlation function is a very "noisy" 
observable: the magnitude of the correlation function is ~ exp(—aLtR), whereas the noise is 
always of order unity. Thus, without any advanced techniques huge statistics is needed. We 
employ here a modification of the multilevel approach presented in Ref. jSJ . 

To summarize our method, let us again consider a lattice of size x Nt x N^, where 
the Polyakov loops are oriented along the t-direction and the correlation function is to be 
measured in the r-direction. One measurement cycle works as follows: 

1. Divide the lattice into Nt sublattices by freezing the following variables: 

(a) gauge link variables Ur and Uz and adjoint scalars x which reside on the (r, 2;)-planes 
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located at t/a = 0,2, ... Nt — 2/ and 

(b) Ut, Uz and x on {t, z)-planes at r = and r/a = Nr/2. 

2. Perform update sweeps on the non-frozen variables within all the sublattices. Be- 
cause the boundaries of the sublattices are fixed, we can generate a valid full lattice 
configuration by choosing any of the configurations from each sublattice, giving a 
total of A^^* lattice configurations. 

3. The Polyakov loop correlation function is measured so that one of the loops is located in 
the interval < ri/a < Nr/2 and the other at Nr/2 < r2/a < Nj.. This gives effectively 
]\fNt measurements; the exponential growth with Nt compensates for the exponential 
decrease in the signal. 

4. In practice, one improved correlation function measurement is (we suppress the z- 
coordinate) 

C(r2 - ri) = P\ri)P{r2) , P{r) = Tr [5(0, r)S{2a, r) . . . S{{Nt - 2)a, r)] , (5.13) 
where S is the average of the 2-link piece of the Polyakov loop within one sublattice: 

^) = 7^ E O^i*^ it + a,T). (5.14) 

1=1 

Clearly, S is an average over A^, P over A^^'/2 and C{r2 — ri) over A^^* measurements. 

The results are further averaged over ri while R = r2 — ri is kept fixed. However, ri 
and r2 should not be too close to the fixed r = and r/a = Nr/2 planes, which would 
diminish the configuration-by-configuration variation and hence the noise reduction. 
Thus, the technique works best when R/a is largest, ~ Nr/2, just where it is needed. 

The parameter A^ is to be optimized for the physical situation in question; the range we used 
is A^ ~ 100 — 300. Naturally, the A^ configurations are far from independent; thus, standard 
update cycles must be performed too. Nevertheless, the method achieves the goal of strongly 
reducing the statistical noise. Clearly, maximizing the number of sublattices achieves the 
highest averaging. However, t-slices of thickness 1 (which would yield A'^^'^* measurements) 
are not very useful, because the variables within blocks would be too strongly locked in place. 

The difference between the method above and the one presented in Ref. |3J is that the 
latter does not implement the freezing of the planes in step 1(b) above. The number of 
measurements is thus cut to A^^*/^; on the other hand a larger optimal A^ tends to compensate 
for this. In our case the existence of the scalar field x general lowers the useful values for 
A^ — indeed, in the broken phase the whole improvement becomes largely unnecessary. 

Let us end this section with another technical comment. Before starting the measurements, 
the system needs to be "thermalised" . That is, one starts from some initial configuration 
''We assume that indices start from zero, i.e., t/a — . . . Nt — 1, for example. 
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Figure 1: Various condensates, as well as the mass parameter y{z), as a function of z/£ at 
a = Ig^ = 2.667, (3 = 9, V = 48^ x 32. In a homogeneous system, the phase transition point 
is at y{z) = tjc ~ 0.15 42 . The classical prediction for the scalar condensate, obtained by 
solving the equations of motion numerically, is also shown. 



(for instance "cold", whereby all fields are frozen to values corresponding to the classical 
minimum of the action, or "hot", whereby all fields have random values) and carries out 
updates until the system has reached a typical "fluctuating" configuration, which should be 
independent of the initial one. It now turns out that with a "hot" start, the thermalisation 
process is anomalously slow. The reason is that the system can contain U(l) vortices and 
anti-vortices penetrating through the domain wall. Since there is no net magnetic flux, all 
the vortices and anti-vortices must eventually evaporate, but microscopic updates are very 
slow in achieving such a global change in practice. With a "cold" start, on the contrary, there 
is no problem with thermalisation, and it is thus the method of choice for this system. 



6. Results 

We now move on to discuss our numerical results. To start with, we show in Fig. ^the actual 
structure of the domain wall for typical parameter values, in terms of local condensates. It 
is seen that, judging by the eye, the width of the domain wall is indeed well determined by 
the input parameter £, as defined through Eq. (|H.3j) . There is quite a difference, though, 
in the magnitudes of the lattice value of the condensate (square boxes) and the tree-level 
value (continuous line). The difference can be explained by the fact that with our choice of 
parameters higher order corrections to condensates are of the same order as the tree-level 
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Figure 2: The static force for various values of a, at /3 = 6, as a function of R/Lr- When not 
visible, the error bars are smaller than the symbol sizes. The volume is V = 24^ for a = oo, 
and V = 32^ X 24 otherwise. The pure gauge value at /? = 6 is taken from Ref. (32) . 



value, so that perturbation theory cannot be trusted. We have checked (in the limit of a thick 
domain wall) that an account of 1-loop effects makes the difference considerably smaller. 

Our basic observable, the force Fp(R) in the central plane z = 0, is shown in Fig. [3 for 
various values of a = ig'^. For a = 0, corresponding to 3d confinining behaviour, the force 
becomes constant (or rather tanh-like) at intermediate distances, representing the 3d string 
tension. Measurements are relatively difficult, since the Polyakov loop correlation function 
decays rapidly with R. As a is increased, the interior of the domain wall starts to look 
more and more like a homogeneous broken symmetry phase, and the plateau moves down, 
eventually disappearing completely. The force then resembles the behaviour observed in the 
3d broken symmetry phase, Eq. (|4.9j) . corresponding to a = oo. 

The quantity we would like to extract from curves of the type in Fig. [3 will be referred 
to as the "string tension". There are, however, two qualitatively different regimes for the 
behaviour observed: for small a, the domain wall is essentially in the symmetric phase and 
we extract the string tension from the plateau in Fp{R) at R <ti Lj.- On the other hand, 
for larger a the domain wall is in the broken symmetry phase, and the string tension is 
extracted from the coefficient of a "linear term" in Fp(R), that is, from Fp(R)/[l — 2R/Lr] 
at i? ~ -Lr/2. It turns out that, as a remnant of the first order transition experienced by the 
homogeneous system at x = 0.20, for certain values of a the central plane can even be in 
one of two metastable branches, exhibiting these two patterns. This situation is illustrated 
in Fig. Etleft). 
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Figure 3: Left: The static force in the central plane in two different runs, at £g'^ = 2, /3 = 6, 
V = 32^ X 24. The system can be in one of two metastable states, resulting in different 
qualitative behaviours for the force. These are referred to as the "thin brane" and the "thick 
brane", because the former has the functional form appearing in the limit i ^ 0, the latter 
in the limit i oo. Right: The "thick brane" static force, divided by (1 — 2R/Lr), at 
different values of L^, compared with results for the homogeneous broken symmetry phase. 
The parameters are P = 6, V = Lr x 24? . 



While a force constant in R at intermediate distances is a signal of non-Abelian 3d con- 
finement, a behaviour linear in 1 — 2R/Lr as R ^ Lr/2 still allows for two different interpre- 
tations: 3d Coulomb phase, characterised by Eq. (|4.9|) . and 2d Coulomb phase, characterised 
by Eq. H4.11() . One can differentiate between the two by approaching R ~ Lr/2, where 



— r = lim — r 

5^ R-.if 9^ 



Fp{R) 



1 - 2R/Lr 



3d U(l) , (6.1) 



16c/2L^ ' 
1 

8l?2.^^ brane 



, 2d U(l) . (6.2) 



Therefore, if we increase the system size L^, the existence of a 2d zero- mode is signalled 
by the force approaching a constant, while for a 3d Coulomb phase the apparent "string 
tension" vanishes in the infinite volume limit as ~ l/L^, and the potential becomes a power 
law, Eq. (|4.7() . These different qualitative behaviours are illustrated in Fig. |3{right). Thus, 
whether dimensional reduction takes place or not can be seen by monitoring the dependence 
of the observable in Eq. H6.1|) on the box size L^. 

Results for the string tension, extracted either from a plateau in the static force Fp{R), or 
as defined in Eq. (|6.1j) . are shown as a function of a in Fig.^ We observe that the behaviour 

®Let us stress again that the existence of a string tension in an Abelian theory is specific to 2d. 
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in turn determines the fit ansatz employed (see the text). 



interpolates between the 3d confinement and 3d Coulomb behaviours. For narrow domain 
walls, or small a<2.0, the system behaves as a 3d confinement phase ("thin brane"). For 
larger a, it is dominated by a 2d zero-mode ("thick brane"), with a string tension scaling as 
~ l/-^brane) in Eq. H6.2() . Note that the system goes over into the "thick brane" regime 
only when the string tension is a factor 3... 4 smaller than in the bulk. Increasing a further, 
the string tension continues to decrease, but the system starts simultaneously to look more 
and more like a 3d Coulomb phase, to which it finally goes over as a — > oo. 

As discussed above, to make a rigorous distinction between a 2d Coulomb and 3d Coulomb 
behaviour, one has to carry out a finite-size scaling study, as in Fig. |3{ right). The corre- 
sponding fj's are shown in the top part of Tabled it is again seen clearly that at a = 2.667 
the string tension is indeed independent of L^, while at a = oo it is not. 

Table ^ also serves to demonstrate that our results are independent of the extents of the 
system in all the directions, for a given lattice spacing. The independence of the lattice 
spacing, for a given physical volume, is demonstrated by Fig. 21 

The discussion so far has been for the force in the central plane, z = 0. In Fig. [5] we show 
how the force depends on z. The pattern is as expected: the behaviour characteristic of the 
properties of the 2d zero-mode is well localised, and outside of the domain wall the dynamics 
is that of the bulk theory. 

All the evidence presented so far supports the conjecture that the low-energy dynamics on 
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Table 1: Finite-volume scaling of cr/g^ as defined in Eq. (|6.1j) . at /? = 6, a = 2.667 (left), and 
a = oo (right). 



the brane is determined by a localised 2d zero-mode. In order to understand the distance 
scale at which corrections to the zero-mode dynamics become important, however, one should 
also determine the mass gap to the higher modes in the central plane. To achieve this, we 
define the quantity 

" - [Fp{R)/(l-2R/L)]' ' ^'-'^ 

where [...]' = d[...]/di?. According to Eqs. 1)4. 11() . 1)4. 13() . the zero-mode does not contribute 
to this observable, and the leading non-trivial contribution at large distances is meff(i?) ~ mi. 

In Fig.iniwe plot a discretised version of 'me{[{R), in the dimensionless form £mes{R). Thus, 
if nil ~ l/-^brane oc l/£, we should obtain a constant value, independent of £. This indeed is 
the behaviour observed, within statistical errors. Moreover, the distance where the plateau 
value is reached is ~ 1, consistent with the scenario that the masses of higher modes still 
also scale as as expected for "Kaluza-Klein type" excitations. 



7. Conclusions and outlook 

The purpose of this paper has been to test with numerical lattice Monte Carlo simulations 
some basic features of a mechanism proposed for gauge field localisation on domain wall 
defects by Dvali and Shifman ^6^. The mechanism relies on non-perturbative dynamics in 
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Figure 5: Left: The static force in various planes, for ig'^ = 2.667, /5 = 9, V = 48^ x 32. 
Right: The corresponding string tensions, ajg^, together with the scalar condensate. For 
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the "bulk" outside of the domain wall, and is thus not easily "proven" to work with analytic 
methods alone. 

Dvali and Shifman considered originally a 3+1 dimensional SU(2) gauge theory, arguing 
that the low-energy dynamics was that of 2+1 dimensional U(l) gauge theory. For technical 
reasons, and since the main physics arguments remain essentially unchanged, we simplified 
in this paper the setting further, and took as the starting point a 2+1 dimensional SU(2) 
gauge theory, with adjoint scalar matter. 

The basic pattern we found can be summarised as follows (cf. Fig. 0]). Suppose the 2+1 
dimensional theory has a large confinement scale A, and consider the dynamics on the brane. 
If the thickness of the brane is of order unity with respect to A~^, the dynamics remains 
the same as in 2+1 dimensions, but the effective string tension decreases rather rapidly with 
the thickness (open symbols in Fig. 0]). As the brane is made thicker, the dynamics becomes 
finally 1+1 dimensional (closed symbols in Fig. 0J), with the "confinement scale" (or string 
tension, ^/a) taking over (almost) smoothly from the 2+1 dimensional value. However the 
transition between the two behaviours actually seems to be discontinuous in our system (cf. 
Fig. Ofleft)). Increasing the width further, ^/a goes down as ^ 1/ V^^hrane, while the mass 
gap to the higher modes goes down as l/Lbrane- An effective 1+1 dimensional description 
with an extremely small ^/a ^ A is only reached for Lbranc ^ A~^, and is then valid only 
at distances 3> .^^brane which implies a kind of a hierarchy problem: there have to be two 
different large scales compared with the dynamical energy scales of the zero-mode system 
(such as masses of bound states, or y/a), namely -Z^brane ^"^^ ^- Moreover, scales such as 
are larger than the masses of some of the higher modes, which means that the contributions 
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from the latter are not suppressed in generic infrared observables. 

How would this pattern change in 3+1 dimensions? If we consider the very scenario 
studied by Dvali and Shifman, then there is in fact a clear difference. The reason is that the 
2+1 dimensional U(l) theory no longer has a string tension associated with a linearly rising 
potential (ignoring the exponentially small value produced by the Polyakov mechanism). 
Thus, repeating our measurements, one would indeed expect to observe a qualitative change 
in the low-energy dynamics, once the width of the domain wall is somewhat wider than the 
inverse of the confinement scale outside the brane: the closed symbols in Fig. 0] would all lie 
on an (almost) straight line, at (almost) vanishing a, signalling a phase transition between 
the bulk-like "thin brane" and localised "thick brane" regimes. 

On the other hand, if we rather consider a non-Abelian case, say a bulk SU(3) theory, such 
that the low-energy dynamics is that of 2+1 dimensional SU(2)xU(l) theory, and probe the 
properties of the non-Abelian part, then the pattern should again be largely similar to what 
we found in this paper. In other words, the domain wall has to be wider than the inverse 
of the bulk confinement scale for a 2+1 dimensional zero- mode to exist, but in general the 
confinement scale of the corresponding low-energy effective theory is smaller than that in 
the bulk only by some numerical factor roughly of order unity, unless the domain wall is 
significantly wider than the inverse of the bulk confinement scale. In the latter case the 
masses of localised bound states (or, say, the string tension ^/a) are expected to scale as 
~ 5^ ~ g'|(£~^)£~^, where gl{jl) is the renormalised 4d gauge coupling. For l/£ much below 
the 4d Yang- Mills confinement scale they are thus of the same order of magnitude or larger 
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than the hghtest of the higher modes, with masses ~ l/£. 

The situation in the 3+1 dimensional case can perhaps be illustrated by noting that, in a 
way, the present mechanism is analogous to the familiar dimensional reduction of 4d Yang- 
Mills theory at a finite temperature T, which occurs only for T above the deconfinement 
phase transition temperature. The role of T is played by The conflict arises because we 
now need a domain wall wider than the bulk confinement scale, corresponding to T below 
the deconfinement phase transition, thus no dimensional reduction. 

It would naturally be very interesting to understand whether the patterns observed in this 
paper change if we go to the physically interesting 4 + n dimensional case, with n > 1. 
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